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Gauge model at finite temperature with massive quarks and at finite 
density on anisotropic lattice 
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Critical properties of QCD and the chiral condensate at finite density are analytically studied on an anisotropic 
lattice in the approximation SU(N) ~ Z(N). Asymptotic behavior of the partition function and its continuum 
limit are discussed. 



In the previous paper Q] we considered the pure 
gluon sector of QCD at finite temperature on lat- 
tice with the anisotropy £ = a a /a T ^ 1 [where 
a a (a T ) is the spatial (temporal) spacing]. At 
present the use of anisotropic lattices £ >> 1 be- 
comes a very popular technique . In this paper 
we include matter fields into the consideration. 
The approximation SU(N) ~ Z(N) used here 
does not cover all features of the SU(N) gauge 
theory, but since the SU(N) and Z(N) theories 
belong to the same universality class , it is com- 
monly believed that the Z(N) degrees of freedom 
are responsible for many important aspects of the 
SU(N) phase structure §. 

1. Z(N) model on anisotropic lattice 

Here we consider lattice QCD in the Hamilto- 
nian limit where all the terms proportional to 1 /£ 
are neglected [|l]||, so the partition function can 
be written as 



Z= cxp{-S G -S F } 

{<x},W},W} 

with the gluon action 



(1) 



where og^ r , v g Z(N), and n T = jpr£- For the 
fcrmion part of the action, we choose the form 
BP] in the approximation £ >> I 



(2) 



r - 7o 



a (x)S 



x'-0 



(3) 



where r is the Wilson parameter (0 < r < 1), rif 
is the number of flavors. 

For gluon fields we use the periodic boundary 



conditions cts,t;, 



and fix the Hamil- 



tonian gauge ctj iT;0 = 1 + S°(Qg - I); fig = 
YI^Iq 1 crg lT; o is the Polyakov loop. Fermion fields 
obey standard antiperiodic boundary conditions. 

Summing over the spatial variables a x - n 
one may easily get 



Sf ff = K T {N T ) J2 tt S nt +n + C.C., 



where the effective coupling k T is defined as 



(4) 



1 - e 



1 + 2e" 



I - e~ 3K - 
I + 2e- 3 "- 



The effective action coincides with the ac- 
tion obtained in B in a strong coupling approx- 
imation where the magnetic part of the action 
is suppressed by factor g~ 2 . On an extremely 
anisotropic lattice (£ — > oo), the similar suppres- 
sion appears due to factor the model may be 
studied beyond the strong coupling area jf],[|,D 
and the limit a T << a a — > may be monitored. 

In the weak coupling region, the effective cou- 
pling k T 

1 I - e- 3N ^ e ~ 3Kr 
Pvt — — in „ , _ 7, — . 

3 1 + 2e- 3Ar - e_3KT 
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limjv.,-— >oo — const leads to the reasonable 
renormalization condition on the bare constant 

We consider the special case r — 1 for the 
fermionic action. Taking into account that 
(S x -p,,x')^ = <W#,x', the action (g) can be rewrit- 
ten as 



with = and 

&x>x = 4 (o-o^)^*'" 1 - (1 + Orm)af ) . 

The chemical potential may be introduced in ac- 
cordance with || as £1 — * e~fl, and, after the 
integration over ipi^ fields and taking into ac- 
count that lim aT ^o(l + a T m) NT = e m / T , we get 
for the Z(3) gauge group 



- s* ft = + — 



V3 



cf> + const (5) 



with .a. = n 

2. / , m , u 

- In ch h ch— 

3 V T T 



1 / ch 3 f+ch 3 f 3 N 
3 n I chf +chf - 4 



0(m, /x) 



1 / chf - ±chf + i^shf ~ 

\ Cll "777 " 7T CIITtt' TT~ Sll 7p 



-7T < 



< 7T. 



In the asymptotic area to » T; fi » T the 
real part of fermion contribution remains essen- 
tial only when to ~ (i but outside this region n 
rapidly disappears. It is also easy to show that 
the imaginary part of fermion contribution disap- 



pears both for ^ — > and (if f >> 1) for ^ 
and 



± 



27T. ^£ m 1. 2£ 

_4_ 7T . JH. _ m. Ai , -|-nri 



±00, 







(G) 



'7 - S eff 



S eff COT " 



The effective action S e 
responds to the extended Potts (N=3) or Ising 

It seems worth noting that the fermion part of the action 
enjoys symmetries r)(m, fi) = r)(m, —/it); <p(n) = — <f>{— /i) 
and rj(m, fi) = m). 



(N=2) model which has been studied in detail in 
[flpfl . Our consideration essentially differs from 
that of [[To) by the specific form of 'external 
sources' r\ and <fi. In particular, the fermion part 
of the action depends on to, fx, T only through the 
combinations of f and ^ . In addition, in our ap- 
proach the parameter to may be arbitrary small. 
All parameters are expressed in physical units, so 
the fermion part of the effective action does not 
change in a T — ► 0; N T — > oo limit. 



2. Mean spin approximation 

Here we apply the mean spin approach pd[ | 
to compute analytically the average value of 
Polyakov loop (f2), which allows us to estimate 
at finite density. Let us introduce fi = 
re 10 = y Y^x ^ x ( mca n spin), where V is the 
volume in lattice units. Adjusting the definition 
of quasiavarages |0| to the considered case 



^qS ( 2Frcos0- 2^ReQ x J 

(<t) V X J 



we may write for the gluon part of the effective 
action 



exp{-S%f) 
with 



_ ■< exp{-5 G } y n 



~< 1 y< 



(8) 



1 2 
— In -< 1 = C~^l k ln l k 



l + 2rcos(6»+2^) 
h = ± U -. (9) 



The gluon part of the effective action up to j 3 
terms can be then written as ^| 

- ^± = 7 r 2 + | 7 V (l + 2r cos 36* - 2r 2 ) . 

V o 



2 It is easy to show that the 'mean spin' method in the 
lowest order in 7 coincides with the simple version of the 
'mean field' method. 
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with 7 = 3 i + 2 e -ik T ■ Therefore we may finally 
write 

S eff 



V 



2inf6r . „ „ S ff 

$ = r/r cos 9 H £=— sin 9 - C — . 

V3 V 



To calculate the partition function we shall 
find the saddle point Oo = roe 1 ® where both 
KcS e ff(9, ro) — minRe<S> e // and |Im<S> e // (9, ro)| 
= min |Im5 e //|. Although at any 9 ^ irn 
ImS e ff [9, ro) 7^ 0, it oscillates in thermodynami- 
cal limit with a very high frequency ~ V sin 9 and 

fv-n ^ m ^ e ff [@i r 0; V) dV' becomes negligible at 
any large V. On the other hand, RcS e ff (9,r ) 
gains the minimum at 77 cos 9 — \rj\, so, tak- 
ing into account that rj > 0, we may put 9 = 
0. Thereby we may conclude that the free en- 
ergy = — ^lnZ does not depend on <fi 
and differs from F(0) simply by 'renormalization' 
77(0) — > r/(/_t). Moreover, the independence F(fi) 
of cj) signals about an implicit charge symmetry 
(D x ^> ft* or 9 <-> —9) in the model which is for- 
mally broken in $. This, in particular, leads to 
<ImO)~™ = 0. 

Therefore, one may write for the free energy 
F ~ -FT max {$(0, r ), $(0,0)}, where the sad- 
dle point ro is defined by a$ g°^ = 0. Failing to 
solve such an equation precisely, we studied it in 
the areas r << 1 and 1 — r << 1 . We find that 
the border line between 'ordered' (ro = (|f2|) 7^ 0) 
and 'disordered' ((|0|) = 0) phases can be roughly 
presented ^as^7 + r/~l. 



= 2 



(^(+)^(+)) + (^(-)^(-) 
l-2e"T(l- 



~)chf 



T dlnZ 
V a dm 
Re(f2) 



1 + e 



3. Conclusions 



This paper considers finite temperature QCD 
on an anisotropic lattice (£ >> 1) in the approx- 
imation SU(N) ~ Z{N). 

We argue that at least in above approximation 
the imaginary part of the action plays a marginal 



3 In fact, the border line with disordered phase has a quite 
complicated structure, however, further details of its de- 
pendance on 7 and r\ may be considered only in a more 
accurate approximation. 



role and the free energy at finite density differs 
from the one at [i = mainly by 'renormalization' 
of the real part of the action: r/(0) — > r}{y). 

The chiral condensate in this model does not 
turn to zero at any S 7^ which indicates chiral 
symmetry breakdown, {tyij}) strongly depends on 
(O) for ^ T — >> 1, so the gluon environment 
plays an essential role in this region. 
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